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Abstract. 

We study the stability of a non-Abelian chromomagnetic vacuum in Yang- 
Mills theory in Euclidean Einstein universe S 1 X S 3 . We assume that the gauge 
group is a simple compact group G containing the group SU (2) as a subgroup and 
consider static covariantly constant gauge fields on 5 3 taking values in the adjoint 
representation of the group G and forming a representation of the group 5(7(2). 
We compute the heat kernel for the Laplacian acting on fields on 5 3 in an arbitrary 
representation of 5(7(2) and use this result to compute the heat kernels for the 
| gluon and the ghost operators and the one-loop effective action. We show that 

_ the only configuration of the covariantly constant Yang-Mills background that is 

stable is the one that contains only spinor (fundamental) representations of the 
group 5(7(2); all other configurations contain negative modes and are unstable. 
For the stable configuration we compute the asymptotics of the effective action, 
the energy density, the entropy and the heat capacity in the limits of low/high 
temperature and small/large volume and show that the energy density has a non- 
trivial minimum at a finite value of the radius of the sphere 5 3 . 



PACS numbers: 04.62. +v, 02.40.Vh, 11.15.Tk 

1. Introduction 

A deep understanding of the physics of quantum Yang-Mills gauge field theory at 
low energies is lacking due to the failure of perturbation theory. At low energies the 
interaction becomes strong which leads to the phenomenologically observed property 
of confinement in quantum chromodynamics. However, the precise nature of the 
physical mechanism responsible for confinement is still not known. From the field- 
theoretical point of view this means that the vacuum of Yang-Mills theory at low 
energies has a far more complicated structure than the perturbative one. One of 
the models of a non-perturbative Yang-Mills vacuum is the chromomagnetic vacuum 
that has been proposed by Savvidy [3D]- Savvidy considered a constant SU{2) gauge 
field (which was necessarily Abelian) in 4-dimensional flat spacetime and showed that 
the trivial zero-field perturbative vacuum is unstable under creation of a constant 
chromomagnetic field. Further, Nielsen and Olesen [27] showed that the Savvidy 
vacuum is itself unstable, meaning that the physical nonperturbative vacuum has an 
even more complicated structure. It has been suggested [2 8) that the real vacuum is 
likely to have a small domain structure with random constant chromomagnetic fields. 
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One way to stabilize the chromomagnetic vacuum is to increase the dimension 
of the space-time. In our papers [5J [B] we studied the chromomagnetic vacuum (for 
general compact gauge groups) in higher dimensions in flat spacetime and showed the 
existence of stable chromomagnetic configurations in dimensions greater than four. 
Another way to stabilize the vacuum is to consider curved manifolds. It is easy to 
see that a positive space curvature can provide an effective mass term for the gauge 
fields on the chromomagnetic vacuum, thus, making the vacuum stable. Of course, 
on curved manifolds the notion of constant chromomagnetic gauge fields has to be 
replaced with covariantly constant fields. 

In our recent paper [llj we considered the finite temperature Yang-Mills theory 
on S 1 x S 1 x S 2 with an Abelian covariantly constant background on S 2 . We showed 
that despite the positive spatial curvature the gluon operator still has negative 
modes for any compact semi-simple gauge group, which means that the vacuum with 
covariantly constant chromomagnetic fields on S 2 does not represent the true vacuum 
of Yang-Mills theory. This happens because on spheres, and, more generally, on 
symmetric spaces, covariantly constant gauge fields cannot have an arbitrary value 
of the chromomagnetic field independent on the spatial curvature; there are severe 
algebraic constraints that force the strength of the chromomagnetic field to be of the 
same order as the spacial curvature. This means that the stabilizing role of the spatial 
curvature is depressed by the destabilizing role of the chromomagnetic field, even on 
curved manifolds with constant positive curvature. 

The present paper is a continuation of the study of Yang- Mills theory on spheres; 
our primary goal here is to extend this analysis to the Einstein universe S 1 x S 3 with 
a non- Abelian covariantly constant background on S 3 with a compact simple gauge 
group G that has the group SU(2) as a subgroup. We refer to the paper [TT] for 
general introduction and notation. 

This paper is organized as follows. In Sect. 2 we briefly describe our model and 
fix notation. In Sect. 3 we derive the one-loop effective action in terms of the heat 
kernels of relevant operators. In Sect. 4 we compute the heat kernel on S 1 . In Sect. 
5 we describe in detail the geometry of S 3 and compute the scalar heat kernel on the 
group SU(2), which is used later for the calculation of the heat kernel for arbitrary 
fields on S 3 in Sect. 6. In Sect. 7 we use these results to compute the trace of the 
Yang-Mills heat kernel (which is the difference between the gluon heat kernel and the 
ghost one) and the effective action. 

The main results of this paper are the calculation of the heat kernel for arbitrary 
fields on S 3 and the proof that the minimal eigenvalue of the gluon operator is positive 
only in a very specific case when the representation of the gauge fields does not contain 
any higher spin , j > 1, representations of the group SU{2) but contains only the spinor 
(fundamental) representation, j — 1/2, of SU(2). In all other cases, that is, when the 
representation of the gauge fields contains at least one representation of SU(2) with 
j > 1, the minimal eigenvalue is negative and the heat kernel grows exponentially at 
infinity leading to the infrared instability of the chromomagnetic vacuum. Of course, 
the gauge fields are taken in the adjoint representation of the gauge group G. The 
question whether the adjoint representation of a compact simple group G, that has 
the group SU{2) as a subgroup, can contain the spinor representation of SU{2), 
is a representation-theoretic problem that we do not solve in this paper. If this is 
impossible then our results indicate the instability of any chromomagnetic background 
in Yang-Mills theory in Einstein universe. 

We assume that this is possible, that is, the adjoint representation of the gauge 
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group G may contain the spinor representation of SU{2), and compute the heat kernel 
and the effective action for this specific case. Of course, it does not make any sense 
to study thermodynamics of an unstable theory. Therefore, in Sect. 8 we study 
thermodynamics of the model for this specific case; we compute the entropy, the heat 
capacity and the pressure. 

The quantum field theory on homogeneous spaces has an extensive bibliography. 
We list only some papers that had an influence on our own approach. First of all, we 
would like to mention the early important papers by Dowker [HI HH EQ1 EH [22] as 
well his papers with coworkers [23] IB ■ A very good source is the excellent review by 
Camporesi [12] and the reference therein as well as the papers [29] [2] [13j [H] [15] [24] . A 
thermal Yang-Mills theory in Einstein universe was also studied in [51] . More recent 
papers [55] [TH] H5] have some overlap with our work since they also studied the heat 
kernel on S* 3 . 

We would like to stress from the very beginning the differences with these papers. 
Our approach is completely different from the papers mentioned above; it is based on 
a fundamental integral representation of the heat semi-group on homogeneous bundles 
as an integral over the isotropy group (see our previous work [4] [71 [8] , for application 
of this approach to quantum gravity see [9]). Also, while most of the authors deal 
with the heat kernel for irreducible (that is, traceless divergence free) tensor fields, we 
compute it for arbitrary fields. 



2. Yang-Mills Theory 

We consider the manifold M = S 1 xS 3 equipped with the standard product metric 
and a compact simple gauge group G. We denote the tensor components with respect 
to the coordinate basis by Greek indices. Let x^ be some local coordinates, A = A^dx 11 
be the Yang-Mills connection one-form taking values in the adjoint representation of 
the Lie algebra of the group G and J- — \F^ W dx^ A dx v — dA + A A A be its curvature 
2-form. Then the classical action of the pure Yang-Mills theory is 

S = [ dvol ^\T\ 2 , (2.1) 
Jm le 

where (ivol = dx g 1 ' 2 is the Riemannian volume element with g = det g^ v , e is the 
Yang-Mills coupling constant and \F\ 2 — g^ a g u P F^J-ap. 

We choose a local orthonormal frame on the tangent bundle TM and denote the 
tensor components with respect to the orthonormal frame by low case Latin indices 
from the beginning of the alphabet. The components of tensors on S 3 are denoted by 
Latin indices from the middle of the alphabet. We introduce the projection tensor h a b 
to the tangent space of S 3 and the Levi-Civita symbol e a bc on S 3 . 

We assume that the background Yang-Mills curvature is parallel V a The = 0. Such 
background on S 1 x S 3 has to satisfy the constraint 

[T abl F cd \=-\ti a [c F^ dh (2.2) 

where a is the radius of S 3 . To satisfy this constraint we assume that the gauge group 
G has the group SU(2) as a subgroup. Then there exist generators Xi taking values 
in the adjoint representation of the Lie algebra of the group G satisfying the algebra 
su(2) 

[Xi, Xj] = —e k ijXk, (2.3) 
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that is, Xi form a representation X of the algebra su(2). Now, the above constraint 
is satisfied by 

Tab = -^e l abXi. (2.4) 
or 

This means that the magnitude of a covariantly constant gauge field cannot be assigned 
arbitrarily but it is determined by the radius of the sphere S 3 , which also means that 
the spatial curvature and the chromomagnetic field are of the same order, I /a 2 . 

In general, this representation is reducible. It decomposes into a sum of irreducible 
representations Xi = diagjY^ 1 © • • • © Y? N }, where ji are the standard non- negative 
integer or half- integer labels of the representations of su(2). The Casimir operator 
X 2 = XiX % is a diagonal matrix 

X 2 = diag (ji + l)I 2n +i © • • • © -3n{3n + l)%v+i} , (2.5) 

where hj+i is the identity matrix of dimension (2j + 1). The labels ji are uniquely 
determined from the adjoint representation of the group G. It is an interesting 
representation-theoretic problem to determine all labels ji given a compact simple 
group G. In particular, we will be interested in the question whether any of the labels 
ji can be equal to 1/2, that is, whether the adjoint representation of the group G 
contains the spinor representation of the group SU{2). We will show below in Sect. 
7 that the stability of the model depends on the answer to this question. We do not 
solve this question in this paper. We simply carry out the calculations for arbitrary 
values of the labels ji. 

We will consider spin-tensor fields taking values in the Lie algebra of the gauge 
group G. We consider a spin-tensor representation E of the spin group Spin(4) with 
generators E & satisfying the algebra 

[£ a6 ,E cd ]=-4<5[ a [c E%. (2.6) 

Recall that Spin(4) = SU{2) x SU{2). That is why the generators 

T 4 = \eijk& k (2.7) 

form a representation T of the subalgebra su(2) 

[T i ,T j ]=-s k ij T k . (2.8) 
Then the matrices 

G t = T t © I x + I T © X, (2.9) 
form the twisted representation T®X of the algebra su{2) and satisfy the same algebra 

[G i ,G j ] = -e k ij G k . (2.10) 

The covariant derivative acting on spin-tensor fields taking values in the Lie 
algebra of the group G is defined as follows. Let T be a general spin-tensor 
representation of the group SO (4), X be a representation of the gauge group G and 
uj ab ^ be the spin connection corresponding to the chosen orthonormal frame. Then 
the covariant derivative is defined by 

V™ x = + l w o6 M T(E o6 ) + X{A,) . (2.11) 
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3. Effective Action 

There exists a gauge such that the one-loop effective action of the Yang-Mills theory 
is given by 16] 

= ~ (logDetii - 21ogDetL ) , (3.1) 

where Det is the functional determinant, the operators Lq and L± have the form 

L = -Ax, Li — —At^x + Q , (3.2) 

Ax = 5 Miy V^fV;f is the Laplacian acting on scalar fields in the representation X of 
SU(2), A Tl8 x = /"Vj lX;f Vj l8X is the Laplacian acting on vector fields in the 
representation X and Q is an endomorphism defined by 

(Q^) a = (R\ - 2F a b ) tp b , (3.3) 

where R a b is the Ricci tensor. 

The most important observation that should be made at this point is that the 
positive curvature acts as a mass (or positive potential) term in the Yang-Mills operator 
L\. While the magnetic field reduces the eigenvalues of the Yang-Mills operator the 
positive Ricci tensor increases them. Roughly speaking, it is the balance of these two 
terms that determines whether or not the Yang-Mills operator is positive. However, as 
we have seen above, the relative value of these two terms is not arbitrary, one cannot 
increase the spatial curvature without increasing the chromomagnetic field at the same 
time. It is precisely this feature that leads to the instability of the model, as shown 
below. 

In order to study the infrared behavior of the system, one has to introduce 
an infrared regularization as in [6]. That is why we introduce a sufficiently large 
mass parameter z so that all operators are positive, which is equivalent to replacing 
the operators L by L + z 2 . The determinants of positive elliptic operators can be 
regularized as follows [31 [TH] . We denote the heat kernel diagonal of the operator L by 
Ui,(t) and its fiber trace Oz,(t) = tr J7t(t). We introduce an arbitrary mass parameter 
A and define the coefficients A k (A) by the expansion as t —> 



Q L (t) ~ (47rf)- 2 e- tA J2 A kW tk ■ (3-4) 

k=0 

Then one can define the renormalized determinant by 

/■OO 1/ 

logDet r0 n(i + z 2 ) = -vol(M) / -e- tz2 &r(t), (3.5) 

Jo t 



where 



Next, let 



er(t) = e L (t)-(«)- 2 e- tt [A^(X) + A^(X)t + A^(X)t 2 } . (3.6) 



&YM(t) =e Ll (;)-2e Lo , (3.7) 

Al M (\) = A L k i(\)-2A L k °(\). (3.8) 



L k V X J "^k 

Then the effective action has the form [TT] 



1 M r dt 



r (1) = --vol(M)j/31og^+ / je- tz 9T«(t)[, (3.9) 
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where = A$ M (0) - z 2 A\ M (0) + ±z 4 A™(0), 

0^ (i) - yM (t) - (4^)- 2 e- tA2 K M (A) + Af M (A)t + Al M {\)t 2 ] . (3.10) 

Finally, there are two simplifications we can make. First, the heat kernel for a 
Laplace type operator L = — A + Q with a covariantly constant potential Q factorizes 

exp(-iL) = cxp(-tQ) exp(tA) . (3.11) 

Second, for the product manifold M — S 1 x S 3 the heat kernel for the Laplacian 
factorizes accordingly 

exp(tA) = cxp(t<9 2 ) cxp(tA S 3), (3.12) 

that is, 

U L {t) = U sl (t)cxp(-tQ)U s3 (t), (3.13) 

where U s (t) and II s (t) are the heat kernel diagonals for the Laplacians on S 1 and 
S 3 . Thus, we obtain finally 

Q YM (t) = U sl (t)tv x [tr Tl exp(-tQ)U% 9X (t) - 2E/| 3 (t)] , (3-14) 

where U§^ x (t) and Uf? (t) are the heat kernel diagonals of the Laplacians acting on 
vectors and scalars on S 3 in the representation X. 

4. Heat Kernel on S 1 

We denote the radius of the circle S 1 by a\. There are two dual representations of the 
heat kernel diagonal: the spectral one 

'"^i^B"')- (4 - 1) 

and the geometric one 

00 / 2 2 \ 

U sl (t) = (4nty^ ]T exp(-^n 2 ). (4.2) 

n— — oo ^ ' 

That is why we will write it in the form 

U s \t) = (4ntr^n^, (4.3) 
where the function Q(t) is defined by 

oo / 2 2 \ 

m = J2 cx p f-^j = 63 {^ 2/t ) > ( 4 -4) 

and 9z{v,q) is the third Jacobi theta function. This function satisfies the Poisson 
duality formula 

t"~ fir 2 



m = yn n {T) (45) 

and has the following asymptotics: as t — > 

f!(t) = l + 2e- 7r2 /* + 0(e- 4 " 2 /*) (4.6) 

and as t — > 00 

= [1 + 2e ~* + ° ( e_4 *)] • (4-7) 
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5. Geometry of S 3 

We consider the sphere S 3 of radius a. However, to simplify notation we set the radius 
o = l. It can be easily reintroduced later from dimensional arguments. Let x % be the 
normal geodesic coordinates on S 3 with the origin at the North pole and ranging over 
(— 7r, 7r). In this section all indices denote tensor components on S 3 . The position 
of the indices on the coordinates will be irrelevant, that 1 S ^ X i — X We introduce 
the radial coordinate r = \x\ = \Jxix 1 , and the angular coordinates (the coordinates 
on S 2 ) 9 l = x l = x l /r, (so that 6i6 l — 1). The round metric on S 3 in geodesic 
coordinates is 

ds 2 = air 2 + sin 2 rdO'de,. (5.1) 
That is, the metric tensor is 

• 2 

sin r 

9ij = Wj + ^2-%' ( 5 - 2 ) 

where 

n tf - Sij - 6,9, . (5.3) 

The sphere S 3 can be viewed as the homogeneous space S 3 — SO(4)/SO(3) or 
S 3 = [SU(2) x SU(2)]/SU(2) = SU{2). The geodesic coordinates on S 3 are exactly 
the canonical coordinates on the group SU (2). Let F{q,p) be the group multiplication 
map in canonical coordinates. This map has a number of important properties. In 
particular, F(0,p) — F(p, 0) = p and F(p, —p) = 0. Another obvious but very useful 
property of the group map is that if q — F(uj,p), then ui = F(q, —p) and p = F(—u>, q). 
Also, there is the associativity property F(oj, F(p, q)) = F(F(uj,p),q) and the inverse 
property F(—p, — oj) = —F(u,p). 

Let Ci be the matrices of the adjoint representation of su(2) defined by 

(C\y k = -2e j ik ; (5.4) 

they satisfy the algebra 

[C i ,C j ] = -2e k ij C k . (5.5) 

Note that the Casimir operator in this normalization is C a C a = —81. Let C = C(x) 
be the matrix defined by 

G = C,x l , (5.6) 

and let R, Y and D be the matrices defined by 

y-i-SLQ, H-i^gj. D.^C. (5,, 

Obviously, RY = YR = 1 and D = Y T R. 

To compute these matrices explicitly we note that 

C 2 = -4r 2 II, (5.8) 

and, therefore, 

C 2n = (-l) n (2r) 2n n, C 2n+1 = {-l) n {2r) 2n C. (5.9) 

Therefore, the eigenvalues of the matrix C are (2ir, — 2ir, 0), and for any analytic 
function 

f(C) = /(0)(J - II) + nl [/(2ir) + f(-2ir)} + ^-C [f(2ir) f(-2ir] , (5.10) 
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tr/(C) =/(0) + /(2*r) + /(-2ir), (5.11) 
det/((7) = /(0)/(2ir)/(-2ir). (5.12) 
This enables one to compute 

Y^-H+^cosrll-i^C, (5.13) 
r 2 r z 

R = 7-n + rcotrn+ -C, (5.14) 

J D = /-n + ncos(2r) + ^^C. (5.15) 
These matrices satisfy the identities 

9ij = S ab Y a t Y b „ gV = 5 ab R\W b . (5.16) 
Therefore, both sets of one-forms 

a\ = Y a l {x)dx\ a a _ = Y i a (x)dx i , (5.17) 

form orthonormal bases on the cotangent bundle of S 3 and the standard metric on S 3 
is bi-invariant. These are the left-invariant and the right-invariant one- forms satisfying 
the identities 

do-% = ±e a bc a b ± A a c ± . (5.18) 

To avoid any confusion we will always use the right- invariant orthonormal basis cr" . 
The right-invariant and the left-invariant vector fields are defined by 

They satisfy the algebra 

[K+,K+]= -2e c ab K+, (5.20) 

[K-,K-]=2e c ab K-, (5.21) 

[K+,K b -}=0, (5.22) 

and form two mutually commuting representations of the group SU(2). The Casimir 
operators of these representations are equal to the scalar Laplacian 

A = S ab K-K- = 5 ab K+K+ . (5.23) 

The left-invariant and the right-invariant vector fields are the Killing vectors of 
the metric generating the whole isometry group of the sphere S 3 . They provide the 
orthonormal bases for the tangent bundle. We will always use the right-invariant 
orthonormal basis K + dual to er" . Note that the geodesic distance between the points 
p and q is just equal to \F(p, q)\. 

The Riemannian volume elements of the metric g is defined as usual 

dvol (x) = g 1/2 (x)dx 1 A dx 2 A dx 3 = sin 2 r dr d6, (5.24) 

where d9 is the volume element on S 2 . The invariance of the volume element means 
that for any fixed p 

dvol(x) = dvo\(-x) = dvo\(F(x,p)) = dvol(F(p,x)). (5.25) 
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We denote the covariant derivative with respect to K+ simply by V . The Levi- 
Civita connection of the bi-invariant metric in the right-invariant basis is defined by 

V a K+ = e c ba K+, (5.26) 

so that the coefficients of the affine connection are 

uj a bc = <j a + (V c K+)=s a bc . (5.27) 

Then 

V a K- = -e b ca D d c K+. (5.28) 
Now, the Riemann curvature tensor is 

R a bcd = -e f cd s a fb , (5.29) 
and the Ricci curvature tensor and the scalar curvature are 

R ab = 2S ab , R = 6. (5.30) 

Let V be the total connection on a vector bundle V realizing the representation G 
of the group SU(2). Then the Yang-Mills connection on this bundle and its curvature 
are 

A = a a + G a , T=^T ab a a + An\ = ^e c ab G c cj a + Ao- b + . (5.31) 

The covariant derivative of a section of the bundle V is then 

V a <p = (K+ + G a )<p, (5.32) 
and the Laplacian takes the form 

A = V a V a = (K+ + G a )(K+ + G a ) = A + 2G a K+ + G 2 , (5.33) 
where G 2 = GiG 1 . Also, the derivatives along the left-invariant vector fields are 

V*- = (5.34) 
where B a = D a b G b . 

We want to rewrite the Laplacian in terms of Casimir operators of some 
representations of the group SU(2). The covariant derivatives V a do not form a 
representation of the algebra SU(2). The operators that do are the covariant Lie 
derivatives. The covariant Lie derivatives along a Killing vector £ of sections of this 
vector bundle are defined by 

Ct = V £ - \a a + {V b ti)e bc a G c . (5.35) 

By denoting the Lie derivatives along the Killing vectors by ICf this gives for the 
right-invariant and the left-invariant bases 

K+ = £ Kt = V a + G a = K+ + 2G a , (5.36) 

JC- = C K - = V K - -B a = K~ . (5.37) 

It is easy to see that these operators form a representation of the isometry algebra 
su{2) x su{2) 

[/C+/C+] = -2e c ab K,t, (5.38) 
[K-,K b ] =2e c ab /C-, (5.39) 
[/C+/C 6 -]=0. (5.40) 
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The Laplacian is given now by the sum of the Casimir operators 

A= 1 -!Cl+ 1 -lCZ_-G\ (5.41) 
where K4=JC±JC± 

We need to compute the action of isometrics on SU{2). Let T a be the generators 
of some representation of the group SU(2) satisfying the algebra su{2) 

[T a ,T b }=-2e c ab T c , (5.42) 
and T(x) — T a x a . First, one can derive a useful commutation formula 

exp[T(x)]T & cxp[-T(x)] = D a b {x)T a . (5.43) 
Next, one can show that 

K+exp[T(x)] = exp[T(x)]T a , (5.44) 

K~ expp»] = T a exp[T(x)]. (5.45) 

Therefore, 

A exp[T(x)] = cxp[T(x)]T 2 = T 2 cxp[T(x)]. (5.46) 
In particular, for the adjoint representation these formulas take the form 

DC b = D a b C a D, K+D = DC ai K~ D = C a D, (5.47) 

A D = -8D. (5.48) 
This immediately leads to further important equations 

exp[^+(p)] expp»] = cxp[T(x)] exp[T(p)] = cxp [T(F(x,p))} , (5.49) 

exp[^- (q)} cxp[T(x)] = cxp[T(< Z )] cxp[T(x)] - cxp [T(F(q, x))} , (5.50) 

where K ± (p) = p a K^. More generally, 

C xp[K+(p) + K~(q)] exp[T(x)} - exp[T(g)] exp[T(x)] exp[T(p)] 

= exp[T(F(q,F(x,p)))]. (5.51) 

Then for a scalar function f(u>) the action of the right-invariant and left-invariant 
vector fields is simply 

exp[K + {p)]f(x) = f(F(x,p)), exp[K-{q)]f(x) = f(F(q,x)), (5.52) 

more generally, 

eMK + (p)+K-(q)]f(x) = f(F(q,F(x,p))). (5.53) 
We introduce the average of the group elements over S 2 by 

Ar(r)= / ^exp[rT(fl)]. (5.54) 
J S 2 4?r 

One can show that this is a group invariant. Therefore, it can be only a function 

of the Casimir operator T 2 . It is determined by the characters of the irreducible 
representations. For an irreducible representation j, this is equal to 

Aj(r) = -i— ]T 008(2^)7,- , (5.55) 
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where the sum goes over integer [i for integer j and over half-integer /i for half-integer 
j. Note that the function Aj(r) is periodic with period tt for integer j and antiperiodic 
for half-integer j, that is, for all integer j: 

Aj(t,r -nn) = Aj(t,r), (5.56) 
and for all half-integer j: 

Aj(t,r-im) = (-l)"Aj(t,r). (5.57) 
We can combine these formulas by writing 

A 3 -(t,r-7rn) = (-l) 2jn Aj(t,r). (5.58) 
For a general reducible representation it is the direct sum of the irreducible ones 

A T (r) = diag{A J - 1 © • • • © A jN }. (5.59) 
Next, we show that for any representation T of SU(2) we have the identity 



cxp(iT 2 ) = (Airb)-^ 2 ^ J dr J d9 rsiarexp f-^j exp[rT(9)], (5.60) 
which can also be written as an integral over R 3 

exp(tT 2 ) = (47rf)" 3/ V / dxrsmrexp(-^—]exp[T(x)}, (5.61) 
where, as usual, r = |x|. First, we compute the following Gaussian integral 

(47rf)~ 3 / 2 e* J dr r sinrexp ^— ^\ 4ir cos(2^r) 

=(H^'" +,, ' + (H 

Therefore, for any irreducible representation j of SU(2) we have 
(47rf)" 3/2 e* y c?r r sinrexp f - ^^ 47rAj(r) 

= ^E{G + M)e-*«" + (i + M) 



e -4M(^-i)t ( 562 ) 



e -4juO-l)t 



= exp[-4j(j + l)t]/,-. (5.63) 

This is nothing but the eq. (|5.60|) for an irreducible representation. The general case 
follows from this trivially. 

We note that the above identity (15.601) can also be written as an integral over the 
group SU{2) 

cxp(tT 2 ) = / dvol(x)U (t,x)exp[T(x)], (5.64) 

JSU(2) 

where 

(r + 27m) 2 



U (t,x)= Y (4 7 ri)- 3 / 2 e t ^ r ^exp 



n— — oo 



41 



(5.65) 



One can show that Uo(t, x) is nothing but the scalar heat kernel on the group SU(2). 
First of all, it satisfies the initial condition 

U (0,x) = 5 s z(x). (5.66) 
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Second, it satisfies the heat equation. This can be shown either by a direct 
computation or, more elegantly, as follows. We have 

= (9 t -T 2 )exp(iT 2 ) 

dvol (a;) (d t - T 2 ) U (t,x)exp[T(x)}. (5.67) 

5(7(2) 

Now, by using the fact that T 2 exp[T(x)] = Aoexp[T(a;)] and by integrating by parts 
we get 

/ dvol(x)exp[T(x)](d t - A )U (t,x) = 0, (5.68) 

JSU(2) 

which gives 

d t U = A C/ . (5.69) 



6. Heat Kernel on S 3 

Our goal is to evaluate the heat kernel diagonal. Since it is constant we can evaluate 
it at any point, say, at the origin. That is why, we will evaluate the heat kernel when 
one point is fixed at the origin, which we will denote simply by U(t;x). By using eq 
(I5.4ip we obviously have 



exp(iA) = cxp (— <G 2 ) exp 
Therefore, the heat kernel is equal to 



(6.1) 



U(t,x) =exp(-tG 2 )* (i,xj , (6.2) 

where ty(t, x) — exp \t (JC 2 ^. + /C 2 .)] S(x) is the heat kernel of the operator (/C 2 . + JC 2 J) . 

First of all, we note that the scalar heat kernel has two important properties: the 
invariance property 

U (t,F(p,q))=U (t,F(q,p)), (6.3) 

and the semigroup property 

dval(p)U (t,F(p,z))Uo(s,p) = U (t + s,z). (6.4) 

S(7(2) 

By using the above method, eq. (|5.64l) . with T = /C+ and T = — /C_ we get 
exp^/C 2 , + tJC 2 _ ) = J dvol (q) dvol (p) U (t, q)U (t,p) 

SU(2)xSU(2) 

xexp[lC + (q)-lC-(p)}. (6.5) 
By using the form of the operators fC^, eqs. (|5.36[) - (|5.37[) . we obtain 

cxp(t/C 2 _ + tJC'i ) = J dvol (q) dvol (p) U (t, q)U (t,p) 

SU(2)xSU(2) 

x cxp[2G(q)] exp[K+{q) - K~ (p)} (6.6) 
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Therefore, by acting on the function \&(s, x) we get 

V(t + s;x) = J dvo\(q) dvol{p)U {t,q)U Q (t,p) 

SU(2)xSU(2) 

x exp[2G(g)]*(s; F(-p, F(x, q))). (6.7) 
We change the variable p by p — F(x, F(q, —z)) so that F(—p, F(x, q)) = z . Then, 

^(t + s;x)= J dvolO) dvol(q)U (t,q)U {t,F(x,F(q,-z))) 

SU(2)xSU(2) 

xexp[2G(q)]V(s;z). (6.8) 
Now, we take the limit s — > and use the fact that ^(0, z) = (^(z) to get 

V(t;x)= [ dvol(q)U (t,q)U (t,F(q,x))exp[2G(q)]. (6.9) 
JSU(2) 

Finally, we compute the diagonal by setting x — 



*(t;0) = / dvol(q)U (t,q)U (t,q)exp[2G(q)]. (6.10) 

JSU(2) 

Now, by using eqs. (|6.2[) . (|5.65l) we get the heat kernel diagonal 

U s \t) = (47rf)- 3 / 2 e* exp [~tG 2 ] S(t), (6.11) 

where 

S(t)= ^2 32tt / dr (r + 27rn)(r + 27rm)A G (r) 



71 ,771 — — OO 



x exp j-^ [(r + 27m) 2 + (r + 2ttto) 2 J j , (6.12) 

and Ag(r) is the average of exp[2rG(#)] over the sphere S 2 introduced above, eq. 
(|5.54p . Note that A G (r) is an even function of r which is periodic with the period 2ir. 
Therefore, the integral can be extended to the interval [— tt,tt]. Thus, the integrand 
is a periodic function with the period 2ir. Therefore, by changing the variables by 
r I— > r — 2ixm and n i— > n + m one summation results in the integration over the whole 
K and we get 

oo 00 

S(t)=2r 3 / 2 ]T j (r + 2™)rA G (r)exp j-^ [(r + 2irn) 2 + r 2 ]X . (6.13) 

— OO 

Next, we change variable by r h > r-\/t — 7m to get 

OO 

S(t) = £ exp L^f) f * ( 2 . 2 - 2^) A G (V* - *») . (6.14) 



n— — oo 



Let us compute the heat kernel diagonal for an irreducible representation j of 
SU(2). By using the explicit form of the function Aj we get 



2j 

oo 



l/*l<j 

OO 

f/r 2 

— r 



^ W = 27TT S M- 3/2eib '°' +1)+1] E (-l) 2j "ex P (-™ 

71 — — OO ^ 

/ 2 2 \ 

2r 2 - 2^- J cos (2p,rVt \ . (6.15) 
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This integral is Gaussian and can be easily computed 
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Uf 



■3/2gt[i0 + l) + l-p 2 ] 



M<3 

x \ 1 - 2^ 2 i - 2- 



00 / 2 2 \ 

exp (-If) 



(6.16) 



This formula is one of main results of this paper. It gives the heat kernel diagonal of 
the Laplacian on S 3 for an arbitrary irreducible representation of the group SU(2). 
At this point it is convenient to introduce a new function 



E (-l) 2j "exp 



^ 2 n 2 



Then 



^ 3(t) = 27TT E (4 J rt)- 3 / a e*yw+ 1 )+ 1 -" a ] {(1 - 2^m 3 {t) 2m' 3 {t)} 



Note that for integer j 

oo 

= = ex p 

n— — co 

and for half-integer j 

oo 

n 3 -(t) = n(t) = ^ (-i) n ex P 



(0,e-^), 



(6.17) 



(6.18) 



(6.19) 



04 (o, e- 7 ^/*) . (6.20) 



By using the duality of the theta-functions we can rewrite these functions in the form 



Q(t) = J-6 2 (0,e- t ) = J- V e" 



(6.21) 



(6.22) 



where the first sum goes over all integers and the second sum goes over all half-integers. 
Here, of course, 8i(v, q) are Jacobi theta functions. 

7. Yang-Mills Heat Trace 

Now we can compute the heat trace for the Yang-Mills theory. By using the results 
for the heat kernel (and restoring the radius a of the sphere S 3 ) we get 



Q YM (t) = (4irt)- 2 n 



t 



e t/a2 ti x \ tr Tl exp 



t 



~(G 2 +a*Q) 



For the vector representation Ti the generator G a has the form 

(G a ) b c = -e b ac ®I x + h b c X a , 



S G 



t 



2 exp ( — — Y 2 Si — 



(7.1) 
(7.2) 
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and the potential term Q has the form 

(a 2 Q) b c = 2h% - 2e b ca X a . (7.3) 
The Casimir operator is then 

(G 2 )\ = -2h b c + h b c X 2 + 2e b ca X a . (7.4) 

Therefore, 

(G 2 + a 2 Q) b c = h b c X 2 . (7.5) 
Thus the above expression simplifies to 

&YM(t) = (4TTt)- 2 n (J^j W (J^j , (7.6) 

where 

W{t) =ehv x a^{-tX 2 ) ^x Tl S G {t)-2S x {t)\. (7.7) 

To compute the function Sg we need to compute the function Ac that is, the 
trace tr^ exp[2G(<7)], where 

2G(q) = C(q) ® I x + Ic <8> 2X(q). (7.8) 
Obviously, these two matrices commute. Therefore, 

tr Tl exp[2G(g)] = exp[2X( g )]tr Tl exp[G(q)]. (7.9) 
The eigenvalues of the matrix C(q) are (0, 0, 2i\q\, —2i\q\)\ therefore, 

tr Tl exp[G(g)] =2 + 2 cos(2|g|), (7.10) 
and we obtain the trace 

tr Tl exp[2G(g)] = 2 [1 + cos (2|g|)] exp[2X(q)]. (7.11) 
Now, by using the form (|6.14l) of the function S(t) we obtain 

n — — rvi \ / v \ 



t 



x 2 cos {^rVt^j tr x exp (-tX 2 ) A x (ry/tj . (7.12) 

Next, by using the explicit form of the function Ax in terms of the irreducible 
representations (|5.55|) we get 

N oo / 2 2 \ 

W(t) = EE E (-l) 2ji "ex P f-^f ) (7.13) 

i=l |/j|<j; n=-oo ^ ' 



x / *L e - r ' [ 2r 2 -2— -- | 2 cus 



(2rVt) cos (2/xn/t) 



— oo 

Now the integrals over r are Gaussian and can be computed; we obtain finally 

N oo / 2 2 \ 

W{t)=Y J E e*C*+i)+i]* 53 (-l)^e X p(-^) (7.14) 

»=1 |/j|<ji n= — oo ^ ' 



1 - 2t(/x + l) 2 - 2—\ e -*(M+i) 2 + A _ 2t( M - l) 2 - 2— ) e-^- 1 
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Also, one can express the function W in terms of the function VLj introduced above 

N 
i=l 

where 

-t( P +i) 2 



Wj(t) = eW+V+i]*l [(1 - 2t(n + l) 2 ) fy(f) - 2tn' J (t)] 

+ [(1 - 2t(ji - l) 2 ) (lj(t) - 2tn' J (t)] e -*(M-i3 2 |. ( 7 . 16 ) 

We will need the asymptotics of the function Wj(t) as t — > and as t — > oo. The 
asymptotics of the function Wj (t) as t — > are 

/22 11 \ 

W 3 (t) ~ (4j + 2) - (8j + 4)i + -j 3 + llj 2 + — i U 2 + . . . . (7.17) 



e -t[(^+i) 2 +(^+i) 2 ] 



By using the dual representation (I6.21j) . (|6.22[) . of the functions flj(t) we can rewrite 
the function Wj in the form convenient for the calculation of the asymptotics as t — > oo. 
After some straightforward but tedious cancellations the result takes the form 

W 3 {t) = JUs/agb-Ci+imn J- [ [(v + If - fa + l) 2 ] 

V \n\<j<\u\ I 

+ [(y - l) 2 - ( M - l) 2 ] e "*[(— 1) 2 +(P-D 2 ]|, (7.18) 

where the summation goes over integer /i and v for integer j and over half-integer fi 
and v for half-integer j (the sum over \x is finite and the sum over v is infinite) . 
Thus, as t — > oo we obtain: for integer j 

Wj(t) ~ Aj2 t 3/2 e -Uf ) ( 719) 

where 

Af n = -j - 1, (7.20) 
and for half-integer j > 3/2 

Wj (t)~j=ff-\)t^e-^r\ (7.21) 



Af n = -i-7- (7-22) 



where 

\j 

Thus, for integer j > 1 

Af n < -2 . (7.23) 
Also, for any half-integer j > 3/2, 

9 
4 



A min<__ ( ?24 ) 
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The only case when the minimal eigenvalue is positive is when j = 1/2; we show 
below that in this case 

19 

K/2 = T - (7-25) 

This is a one of the main results of this paper. It tells us that the minimal 
eigenvalue is positive only in the case when the representation X does not contain 
any higher spin representations of SU(2) with j > 1 but contains only the spinor 
(fundamental) representation of £{7(2) with j = 1/2. In all other cases, that is, when 
the representation X contains at least one representation with j > 1, the minimal 
eigenvalue is negative and the heat kernel grows exponentially at infinity, at least as 
e 2t , leading to the infrared instability of the chromomagnetic vacuum. 

Now, we can compute the effective action. First, by using the asymptotics of the 
function Wj (t) as t — > we get the asymptotics of the heat trace 

+2 



Q YM (t) ~ (47rf)- 2 | Co + Ci^ + C 2 ^ + ■ ■ •} , (7.26) 



where 



N 



C =^(4ji + 2), (7.27) 

i=i 

N 

Ci = £(-8* - 4), (7.28) 



i=l 
N 

a 



E (f^ 3 + n 3i + yJi) ■ ( 7 - 29 ) 

Therefore, (3 = (47r)~ 2 [C 2 — z 2 a 2 C\ + ^z 4 a 4 Co] /a 4 . Now, the renormalized heat 
trace becomes 



i ., / t \ „ r ( t \ „ ft 



e^ M (t) = (iTTty 2 ^n^-jw^j- r ym j , (7.30) 

where 

Rym (t) = e- ty? j<7 + (d + C A 2 ) i + (^C 2 + d\ 2 + \c \^j t 2 } . (7.31) 



Thus, the one-loop effective action is 



r - - n 



C 2 - z 2 a 2 C\ + -z 4 a 4 Cv 



log^ I •!•<:.<•;. (7.32) 



where x = a/ai and (after rescaling the integration variable t — > to 2 ) 

f°° dt 2 2 

$(a;) ^ e~ ta 2 {fi (x 2 t) T4^ (t) - R YM (t)} . (7.33) 



/o 

The total effective action (including the classical action) is then 



7T f 87r 2 cr 
8x\ e 2 



C 2 - zVd + iz 4 a 4 C 



log§f 'l" : - r: ' (7.34) 



where a = —trX 2 = X)i=i + This formula is another important result of this 
paper. It gives an exact integral representation of the infra-red regularized one-loop 
effective action for pure Yang-Mills theory in Einstein universe. 
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Recall that z is an infrared regularization parameter introduced to ensure 
convergence of the integral (|7.33l) at infinity. Eventually we need to take off the 
infrared regularization, that is, to take the limit z 0. The convergence of this 
integral in this limit depends on the asymptotic properties of the function W(t) at 
infinity. By using the asymptotics of the function W(t) at infinity obtained above we 
see that as a function of the infrared regularization parameter z the effective action 
T(z) is analytic for Rez 2 > — A m i n /a 2 , where 

A min = min A™ in , (7.35) 

Ki<N 



and has a branch singularity at z — \J — A m j n /a. Therefore, if A m i n > the effective 
action has a well defined regular value as z — > 0, but if A m ; n < the effective action is 
singular in this limit and the model is unstable at low energies. 

Recall that the stable configuration only occurs when the adjoint representation 
of the gauge group G contains only the spinor representation of SU(2). In principle, 
for any given group G containing the group SU(2) as a subgroup we should be able 
to determine all the labels ji of the representation X, which might be an interesting 
representation-theoretic problem. If one could show that among the labels ji there 
cannot be any label j = 1/2 then, we would simply conclude that the model is unstable 
for any compact simple group, a very strong assertion. 

However, we do not do this in this paper. We will simply assume that this 
is possible and proceed with the calculation. That is why we will consider further 
the case when the representation X contains only spinor representations, that is, all 
ji = 1/2. In this case the function W simplifies 



W(t) = N J] (-l)"^ 2 " 2 /* 

n— — oc 

= iv| (2 - 9t) Q(t) - 4ttt'(t) 
Therefore, as t — > 



2-9^-4 



TT 2 n 2 



-t/2 



a -t/2 



(2 - t) O(t) - UQ!{t) 



3*3/2 1 



11 



W{t) - N 4 - 8t + — i 



so that 



C Q =AN. 



Ci = -8JV, 



C 2 



11 



-TV. 



(7.36) 



(7.37) 



(7.38) 



Further, by using the dual representation of the function Q we obtain 



W(t) = N^=t 3 ' 2 e~ t{n+1 ' 2 ^ \ (n - l)(n + 2)e~ t ' 2 + n(n + l)e 

* n— 

Therefore, as t — > oo 
W(t) 



t3/2 



48 



j Vt 3/2 e -19t/4_ 



(7.39) 



(7.40) 



It is convenient to introduce the function 
1 



W(t) 



Then as t -> 



W(t) 



lit 



101 2 
~^t 2 



(7.41) 



(7.42) 
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and as t — >• oo 

AO 

W{t) ~ -^Lt 3 / 2 + ■■■. (7.43) 

In this case the effective action is well defined even in the limit when the infrared 
regularization parameter z is set to zero. Therefore, by taking the limit z — ¥ and 
choosing the parameter A by A = \/l9/(2a) we obtain 

r -5»{£-T*(s«V)-*4 

where 



i(a,) = * e" t19 / 4 jfi (x 2 i) # (t) - 4 - lit - ™r 
Note that as x — > this function approaches a well defined constant 



(7.45) 



$o=/ ^ e - tl9 ^\w(t)-A-nt-^-t 2 \. (7.46) 



* 3 

Therefore, we can split the integral in two parts 

$(aO = *o + *iO»0. ( 7 - 47 ) 
where (after rescaling i t/x) 

* x {x) =x 2 f o °°§ exp [Si (xt) -1}W (£) . (7.48) 

By using the asymptotics of the function [Sl(t) — 1] it is not difficult to see that 
there is a critical value 

\/l9 

x c = ^vr. (7.49) 

If x < x c then the integral &i(x) is exponentially small. Indeed, in the limit x — >• 
the function [Sl(xt) — 1] is determined by the first exponential term, (|4.6[) ; therefore, 



7 f°° dt [ 1 /19 7T 2 \1 - /t\ 



(7.50) 



The main contribution to this integral comes from the neighborhood of the point 
to = 27r/vl9. Therefore, by using the asymptotics of the function W(t) (|7.43|) we 
obtain 

* / n 96 5/2 ( VT9n\ 

®i(x) x 5/2 exp . (7.51) 

7T \ X I 

By rewriting the integral (|7.45j) in the form 

* W _ jT * {a ( t , (±) - 4 - n± - Mil} , (7 . 52) 

one can obtain the asymptotics of the function as a; — > oo 

$(x) - ra 4 , (7.53) 
where v is a positive constant defined by 
f°° dt 

u = ij ^[0 (*)-!]. (7.54) 
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It can be computed exactly v = JrC(4) = ^q, with ((s) the Riemann zeta function. 
Thus, 

$(x) ~ ^x 4 . (7.55) 

By expanding the function W(t) in a power series one can also obtain the expansion 
of the function $(x) in inverse powers of x. 



8. Thermodynamics of Yang-Mills Theory 

In this section we investigate the entropy and the heat capacity of the model for 
the stable configuration of the background fields containing only one representation 
with j = 1/2. The temperature is related to the radius of the circle S 1 by 
T = l/(27rai) = x/{2ira) . The volume of the space is the volume of the sphere 5 3 , 
V = 2ir 2 a 3 . For a canonical statistical ensemble with fixed T and V the free energy 
F = E — TS is a function of T and V defined by the total effective action F = TV , 
with T = T(o) + Then the entropy, the heat capacity at constant volume and the 
pressure are determined by the derivatives of the free energy 

c, 9(TT) d 2 (TT) d(TT) 

S-— C v --T^^, (8 .i) 

It is easy to see that neither the classical part T( ) nor the renormalization logarithmic 
term in the combination TT depend on the temperature. Therefore, the entropy 
and the heat capacity do not depend neither on the classical term nor on the 
renormalization parameter /x. Thus, the entropy and the heat capacity (per unit 
volume) are given by the derivatives of the function $, 

By using the asymptotics of the function we can compute now the entropy 
and the heat capacity of the gluon gas: at low temperature and small volume, as 
aT — > 0, we have 



38 Ar a -5/2 T i/2 exp [_v^| (83) 




And at high temperature and large volume, as aT — > oo, we obtain 

8 ~ A 7r 2 7VT 3 ; Cy „ JL^NT 3 . (8.5) 
45 15 

One can also compute the energy density e — E/V = T(T + S)/V to get 
N f 6tt 2 11 /4 



327r 2 a 4 

and the pressure 



%- ~ T log (la" 2 " 2 ) + xi ' (x) ~ ' (8 ' 6) 
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Notice that both the pressure and the energy density depend on the renormalization 
parameter fi. Further, this gives the equation of state of gluon gas 

1 UN , , 



-^-log -«V -$o ■ (8.9) 



3 967r 2 a 4 ' v ' ' 

In the limit of low temperature and small volume as aT — > the energy density 
has the form 

N (6tt 2 U ( 4 
6 ~ 32^1 \— ~ Y og ^19' 

In the limit of high temperature and large volume, as aT — ¥ oo, the energy density is 

•-aMS-T'^)}^- <8 ' 10 > 

Of course, in this limit the gluon gas behaves like the "colored" photon gas that has 
2N as many degrees of freedom. Recall that for the photon gas P = |e = j^tt 2 T 4 
and s = ±Tr 2 T 3 . 

45 

Notice that at a fixed temperature T as a — > the energy density e — > +oo and 
as a — > oo it goes to a positive constant. Also, both functions ()8.9[) and (|8.10[) have 
minimums at some different values of the radius a, 

I /12tt 
-6i j2 exp — - 
(i \U e 

where 



ai,2 = -&i,2cxp ( — — ) , (8-11) 



/1 9 / l\ /19 / l \ 

jexPy-Z)' fe 2 = y T exp(---$ J. (8.12) 

Therefore, it is reasonable that the energy density has a minimum at a finite non- 
trivial value of the radius ao(T) — ~b(T) exp ^jf fr^ between a% and 02, which can be 
determined numerically. Obviously, this is a purely non-perturbative effect; it indicates 
the existence of a non-trivial vacuum with the magnetic field and the curvature of the 
order ~ fj? exp ^— jf 



9. Conclusion 

The primary goal of this paper was to study of the low-energy structure of the Yang- 
Mills vacuum. We assumed the chromomagnetic nature of the vacuum with covariantly 
constant chromomagnetic fields, which was well known to be unstable in flat space. We 
noticed that the potential term of the gluon operator L\ = —At^x +Q has the form 
Q a b = R a b — 2F a b- Therefore, a large Ricci curvature increases the minimal eigenvalue 
and a large magnetic field decreases it. Therefore, to make the gluon operator positive 
one needs a large Ricci tensor and a small magnetic field. Moreover, one needs the 
ability to independently change the magnitudes of the curvature and the magnetic field 
to make this work. 

Of course, to be able to carry out the calculations one needs some degree of 
symmetry. So, we assumed that the chromomagnetic field is covariantly constant and 
the spatial curvature is constant and positive. In four dimensions there are not many 
choices of spaces of non-negative constant curvature, only products of spheres and 
circles; the non-flat ones are S 1 x S 1 x S 2 and S 1 x S 3 . The case of S 1 x S 1 x S 2 was 
studied in our previous paper [IT] and in this paper we studied the case of S 1 x S 3 . 



Thermal Yang-Mills Theory In the Einstein Universe 



22 



One of our results is the calculation of the heat kernel of the Laplacian acting on 
arbitrary fields on S 3 in arbitrary representation of SU(2). By using this result one can 
compute the one-loop effective action of any field-theoretic model. We apply it to the 
calculation of the effective action of Yang-Mills theory on Euclidean Einstein universe, 
S 1 x S 3 . Our main result is the proof that generically the gluon operator almost 
always has negative eigenvalues. The only case in which the gluon operator does not 
have negative eigenvalues occurs when the fundamental (spinor) representation of the 
group SU (2) can be embedded in the adjoint representation of the gauge group (that 
contains the group SU(2) as a subgroup). If one can show that this is impossible, 
then our results prove the instability of the chromomagnetic vacuum of Yang-Mills 
theory in Einstein universe for any compact simple gauge group. We intend to study 
this representation-theoretic problem in a future work. 

We confirm the conslusion of our previous work [llj that to stabilize the 
chromomagnetic vacuum at lower energies one should consider non-constant magnetic 
fields on non-compact spaces. Covariantly constant magnetic fields on compact 
symmetric spaces are too rigid, they are completely determined by the spin connection 
and are of the same order as the spatial curvature. This makes it impossible for the 
gluon operator to be strictly positive. 
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